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射流 气泡 的 被 动 声 源 成 像 方 法 
初 宁 刘 竞 好 厅 世 杰 吴 大 转 
(浙江 大 学 能 源 工程 学 院 化 工 机 械 所 ， 杭 州 310027) 


摘 要 被动 声 源 成 像 方法 近年 来 越 来 越 成 为 研究 气泡 声 源 的 有 效 手段 之 一 . 本 文 所 阐述 的 气泡 声 源 成 像 方 法 ,不 仅 能 
从 有 限 数量 的 水 听 器 测量 信号 中 对 气泡 声 源 进行 高 分 辨 率 的 空间 定位 ， 还 能 重 构 出 气泡 声 源 的 平均 功率 分 布 情况 . 为 此 ， 
本 文 首先 需要 建立 一 个 合适 的 气泡 声 源 模 型 、 进 一 步 简 化 气泡 声 源 在 空间 传播 的 多 径 效 应 ,以 便 充 分 利用 阵列 信号 处 理 的 
优势 ; 然后 , 本 文 简要 介绍 基于 波束 形成 、 最 小 方差 无 失真 响应 和 子 空 间 分 解 等 信号 处 理 技术 的 气泡 声 源 成 像 方法 ; 并 和 针 
对 研究 的 射流 气泡 特点 、 流 速 和 喷 管 性 状 等 条 件 ， 还 讨论 了 有 限 数量 的 水 听 器 阵列 对 声 源 成 像 结果 的 影响 . 最 后 ， 本 文通 
过 仿真 模拟 说 明了 被 动 声 源 成 像 方法 可 以 揭示 气泡 声 源 在 气泡 形成 、 生 长 和 脱离 喷 管 等 阶段 的 声 源 平均 功率 高 分 辩 率 空 
间 分 布 情况 . 该 方法 能 够 有 效 “ 和 看 见 ” 气泡 声 源 、 弥 补 了 高 超速 相机 拍摄 气泡 图 像 的 不 足 . 
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(Institute of Process Equipment, College of 到 人 eering, Zhejiang University, Hangzhou 310027, China) 


Abstract Nowadays passive acoustic wgi g has become one of the most effective methods in 
studying the sound emission of jet bu St low airflow rates. This paper presents a proof of concept 
that our proposed approach base Sx atial filtering techniques can not only localize bubble acoustics, 
but also reconstruct acoustic by only using the limited measurements at hydrophone array. First 
of all, an appropriate ac odel of jet bubbles should be investigated carefully, and the acoustic 
propagation should I appropriately in order to fit for arraysignal processing. Then the 
state-of-the-art spatial filters such as Beamforming, Capon and Music M8thods are extended for bubble 
acoustic imaging. And the deployment of hydrophone array is derssed according to the conditions 
of jet bubbles, airflow rates and sparged tank. Finally acoustiNitaging results from simulations can 
show vividly acoustic power distributions with high plation and moderate dynamics during 
the jet bubble formation, detachment and os 
e 


on this proof of concept, limited number of 
hydrophone array can be effectively used in undewwd 


experiments. Moreover our proposed approach 
allows to visualize the bubble sound CR with respect to bubble images from high speed 


camera. SO 


Key words jet bubble acoustics; passive acoustic imaging; source localization; array signal process- 
ing 
0 引 富 


被 动 声 源 成 像 是 一 种 利用 非 接触 式 传感器 阵列 。 理 方法 不 同 ,被 动 声 源 成 像 通 常 能 够 分 为 近 场 声 成 


来 获得 声 源 定位 和 声 强度 绘制 的 先进 技术 .被 动 声 
源 成 像 能 够 对 振动 和 转动 、 摩 擦 和 爆破 、 喷 气 和 射 
流 等 过 程 中 产生 的 声 源 提供 深入 的 机 理 分 析 和 特 
性 统计 0-。 近 年 来 ， 高 分 辨 率 的 近 场 被 动 声 源 
成 像 方 法 已 成 为 水 下 低速 射流 气泡 声 源 研究 的 突破 
性 技术 之 一 -9 根 据 声 源 的 物理 机 理 和 求解 的 数 


收 稿 日 期 :2017-01-13; 修订 日 期 ，2017-04-19 


像 方法 (NAH) 即 所 谓 的 “ 声 相机 ?n9， 基 于 “ 波 
东 形 成 ”的 空间 滤波 器 方法 叫 ， 以 及 基于 “ 解 卷 
积 ” 的 反问 题 推断 方法 - 本 .考虑 到 成 像 性 能 能 
稳健 性 、 计 算 的 快速 性 和 水 听 天 阵 列 布置 的 简易 
性 ， 本 文 重点 关注 基于 “波束 形成 ”被 定 声 源 成 像 
方法 。 
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(a) 气泡 生源 和 传感器 


1 气泡 声 源 的 传播 波 模型 
在 建立 模型 之 前 ， 必 要 的 假设 如 1) 气 
泡 声 源 假设 是 非 相关 的 单 极 子 ， 虽 嫁 这 是 一 种 理想 


的 假设 情况 ,但 是 从 统计 观点 i 头 量 的 单 极 子 
的 整体 性 能 能 够 在 某 种 程度 出 实际 声 源 的 特 


性 . 2) 水 听 器 假设 是 全 向 的 同一 增益 的 . 3) 环境 
背景 噪声 假设 是 加 性 高 斯 白 噪声 . 4) 声 源 传播 服从 
球面 波 、 多 径 效应 主要 考虑 一 阶 的 反射 和 折射 . 

在 图 1(a) 中 ， 气 泡 声 源 和 水 听 器 阵列 布 署 在 
同一 个 平面 内 。 设 有 天 个 非 相关 性 的 单 极 子 信 源 
一 国有 三 ] 。 这 天 个 信 源 位 于 未 知 位 置 
P* 于， 其 中 态 代 表 声 源 大 的 
和 村, [代表 并 轩 和 于 . 为 了 重 构 这 大 个 过 
要 使 用 数量 大 于 K 的 M 个 水 听 器 (M > K) 位 于 
已 知 的 位 置 P= [ 丫 ,… 成 ,,… 友 4] 。 为 了 提高 空间 
分 辩 率 , 信 源 平面 被 离散 化 成 为 N 个 足够 小 的 网 格 
点 s = [s1,…sn,…sN]!'， 且 NN > KK， 其 网 格 点 的 
位 置 已 知 为 已 = [p1,…pn,…pN] ,而 每 个 网 格 点 
可 以 被 看 作 是 潜在 的 声 源 单 极 子 。 所 以 s* 和 s,p* 
和 了 的 关系 可 表示 为 ; 


Sk = Sn 
px — pn 


因 s 的 信号 长 度 N 远大 s 的 非 零 分 量 个 数 KK, 所 以 
5 一 [0 …0,…， , 5105 0， 
站 工 


sn #0 (1) 


(2) 


© ne 
图 1 气泡 声 源 测量 与 < 
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Fig. 1 Acoustic bubble sources ~ 
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rray signal processing 


因此 , s 是 一 个 天 维度 的 N 长 稀 梳 信号 , 即 s 包含 
了 较 少 的 天 个 非 0 分 量 和 较 多 的 N - 天 个 0 分 量 . 
因此 , 重 构 s* 等 价 于 重 构 K 维度 稀 朴 信号 s。 
在 建立 了 声 源 模型 之 后 , 信号 处 理 过 程 如 图 1(b) 
所 示 . 设 第 m 个 水 听 器 测量 到 7 个 采样 时 间 信 号 
点 ,并 将 其 等 分 为 了 个 时 域 信号 段 , 每 彼 有 工 个 时 
间 信 号 点 ,满足 Er Tx 工 . 设 zmlt) 为 第 m 个 水 
人 网 最 (; [0，… ,也 ) 的 第 + 个 时 间 点 
,i]) 所 接收 信号 向 量 . 由 于 时 间 
RN “ zim(t) 通常 具有 宽频 带 信号 特性 , 所 以 在 
只 从 间 段 内 通过 离散 傅立叶 变换 把 xm(b) 转化 到 
频 j nn 并 宽带 频谱 独立 分 为 工 个 子 频带 
(es [如 ). 那么 对 于 所 有 的 M 个 水 听 器 , (有 
代表 全 部 的 上 总 搂 收 信号. 由 于 个 频带 信 呈 可 以 
独立 并 行 的 处 理 ， 后 文中 省 略 有 ,因此 ,zi( 所 ) 能 够 
表示 为 : 


%i= A(P)si+ ei (3) 


其 中 ,4(P) = [a(p1),…a(pn),… ,a(pN)] 代表 水 听 
器 阵列 对 于 声 源 s 的 指向 矩阵 ; a(pn) 代表 水 听 器 阵 
列 对 于 位 于 pn 位 置 的 声 源 sn 的 指向 向 基 ; 对 于 球 
面 传播 的 声 源 ，a(zn) 表达 为 : 


= exp (—j27 fiTn,1) . 


Tn,1 


| Tr (4 
exp (—j27 fiTmn,m) exp (—j27 fiTn, Mm) 
Tn,m Tn,M 
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其 中 , rnm 代表 从 声 源 到 水 听 融 m 的 距离 ,Thm 
代表 历经 rm 所 用 的 时 间 ; 后 文中 ，a(pn) 被 简化 
an。 特别 的 ，rnm 应 该 考虑 到 反射 和 折射 等 多 径 效 
应 , 在 作者 文章 中 讨论 过 Da， 

综 上 所 述 , 公式 (3) 中 的 声 源 传播 模型 是 一 个 
从 M 次 测量 中 求解 NN 个 未 知 数 的 从 定性 线性 方程 
(N > M)。 因此 通常 情况 下 ， 需 要 增加 水 听 融 个 数 
来 降低 求解 的 不 确定 性 . 但 是 受到 实验 条 件 限 制 , 水 
听 融 数量 不 可 能 明显 增加 . 并 且 , 被 动 声 源 定位 的 目 
的 是 定位 声 源 位 置 和 绘制 声 功率 分 布 , 那么 公式 (3) 
可 以 通过 波束 形成 技术 [1 转化 成 为 


yn = i Riin (5) 
其 中 , yn 是 波束 形成 方法 所 估计 出 来 的 在 位 置 pn 处 


的 声 源 功率 及- 是 测量 声 源 的 互 功率 谱 矩 阵 ( 协 广 
关 短 隐 ), 而 如 Te 是 流 来 形成 对 于 声 源 sm 的 


指导 向 量 ,. (-) 是 共 斩 转 置 算 子 . 由 于 声 源 s J 


非 相关 的 单 极 子 ， 并 且 与 表 人 
的 , 因此 丸 . 可 以 进一步 推导 出 A 
> 


= AR,A’+ A 
黎 ' (6) 
> TnQana a 


其 中 ， 当 时 间 分 段 数目 足够 大 (7 > 1), Bls2] 
工 

-2 能 够 快捷 的 近似 计算 ， 且 E[] 是 数学 期 望 
算 子 .Rs。 = Blss'] 待 测 声 源 的 互 功率 谱 ， 由 于 声 
源 s 的 元 素 是 相互 非 相关 的 ， 因 此 R。 = Blss] = 
diag [xz1,.… ,Xn,.…ZN]， 其 中 zn = 五 [sus/， rR 
和 pn 的 声 源 平 均 功率 diae[] 是 到 和 了 对 

算 子 ， Dy 是 M x M 维 的 单位 盾 阵 ; 于 是 将 公式 

带 入 到 公式 (3) 得 到 ， 


N 
yn = (> me Gn + 0 Gn (7) 
y= 


Kh 波束 形成 指向 向 最 fr = 和， 并 且 tn = 
es 所 以 公式 (7) 得 到 
N 
lol ,1 
ra 


最 后 , 对 于 NN 个 洪 在 的 声 源 (网 格 ), 我 们 用 算 
阵 形式 改写 公式 (8), 从 而 得 到 声 源 平 均 功 率 的 传播 
模型 为 


y= Or+ola (9) 


其 中 , z = [zi1,…zn,…zN] ”代表 未 知 的 声 源 功率 
向 量 ; 根据 声 源 功率 的 定义 ， 以 及 无 关 性 单 极 子 信 
源 的 假设 ,功率 x 继承 了 信号 s 的 K 维 稀 朴 性 信 
号 的 特点 . 因此 狭义 上 的 声 源 成 像 能 够 把 从 公式 (3) 
中 重 构 声 源 s 转化 为 从 公式 (9) 中 重 构 声 源 平均 功 


率 zo2 = Eleied 代表 背景 噪声 的 功率 (方差). 向 
1 1 1 
量 4 = | 代表 水 听 器 


阵列 在 不 园 网 格 位 置 上 的 衰减 因子 . 矩阵 C 代表 水 
听 器 阵 乔 允 声 源 功率 的 指向 矩阵 ， 甚 元素 co < C 


Ky 


1 1 
ep (10) 


f 2 
. (ram 一 rnm) } 
0 


其 中 , ay 和 @ 症 传感器 阵列 在 网 格 位 置 ps 和 p, 上 
的 指向 向 前 ,其 定义 在 公式 (4) 中 ,ri 是 从 网 格 n 
到 敢 折 器 m 的 距离 ,i 代表 在 第 1 子 频带 上 的 频率 
全 MM 是 水 听 器 的 总 数量 ; 公式 (10) 揭示 出 功率 指 
向 矩阵 C 的 重要 性 质 ， 


_ | >0 for gn 
| 人 (11) 


事实 上 , co 可 以 被 解释 为 邻 域 网 格 9 对 计算 网 
格 n Bn 例如 公式 (8) 中 


= emt Te ,| 意味 着 , 在 网 格 n 处 的 波束 


形成 功率 是 N 个 潜在 声 源 功 率 的 加 权 平均 入， 该 权 
重 can 由 阵列 的 结构 和 阵列 与 声 源 之 间 的 位 置 决 定 ; 
在 理想 条 件 下 ，cqn ~ bqn 接近 Dirac 函数 、C ~ Tn 
退化 为 单位 矩阵 , 从 而 公式 (6) 简化 为 y = z+o?1， 
从 中 非常 容易 的 求解 出 声 源 功率 z。 这 从 侧面 说 明 ， 
合适 的 阵列 结构 能 够 大 大 有 助 于 获得 高 分 辨 率 和 准 
确 的 声 源 成 像 结 果 . 
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(a) 单 极 子 声 源 在 线形 均匀 分 布 阵列 的 信和 号 接收 (b) 空间 到 达 角 的 定位 能 力 比 较 


图 2 极 子 声 源 在 均匀 线 阵列 的 分 辨 率 比 较 


Fig. 2 Comparison of spatial resolution for monopole source signal at linear uniform sensor 


综 上 所 述 , 相对 于 声 源 传播 公式 (3), 声 源 功 率 ”其 中 , mo 代表 声 源 sn 0 GC 2(a) 


三 。 
~y rz 是 K2 i 量 。 攻 1 
I 果 是 把 y 当 作 z 的 估计 值 ， 虽然 计算 简单 和 稳健 ， Ry ie I 

站 但 是 空间 分 辩 率 往往 比较 低 , 为 此 ， 基 于 最 小 广 | , 

Tap . [xA sing 

己 无 失真 响应 1 (Capon 滤波 技术 ) 的 声 源 成 像 2 全、 nD) 

之 ”能够 有 效 提高 分 辩 率 并且 合理 的 估计 出 点 YnS ~ Ta | (15) 
CO 值 , Capon 滤波 的 主要 思想 是 是 保持 蘑 -二 入 和 的 "|| 

已 伟 感 居 隆 列 的 测量 并 且 和 多 甬 的 测 


该 式 中 含有 Dirichlet 函数 , 其 角 分 辨 率 可 以 表示 为 : 


功率 , 因此 Capon 对 声 源 成 人 的 和 全 为 


【ew | Capon _ 

A 让 a 和 on . Ac Hp 和 
PP 进一步 ， 为 了 得 到 空间 上 的 超 高 分 辨 率 ， 基 于 

> 子 空间 分 解 (MUSIC) 的 声 源 成 像 方法 能 够 从 测 as 最 小 方差 无 失真 响 
人 量 信 号 的 互 功 率 谱 中 通过 奇异 值 分 解 出 相互 正 交 的 人 。 该 仿真 
.二 ” 声 源 空间 Bs 和 噪声 空间 及， 当 位 于 网 格 n 的 声 源 独立 的 高 斯 白 噪 声 的 声 源 (K = 2) 从 位 于 


- 空间 与 噪声 空间 的 正 交 性 越 大 时 , 该 网 格 n 处 越 有 ~ 和 +5? 的 两 个 方向 传播 的 声 信号 ， 且 振动 频率 
只 可 能 存在 声 源 信号 . 因此 WA 为 200 Hz( 低 频 声 源 )， 平 均 功率 均 为 单位 1; 信号 比 

率 为 : 为 0 dB, 即 声 源 功 率 和 背景 噪声 功率 相当 ; T=10000 
2 (13) ”个 时 间 采 样 点 ; 声 源 平面 被 等 分 成 为 V=361 个 网 格 ， 


2 
ah.BEell 


2 射流 气泡 声 源 成 像 仿真 

在 图 2(a) 中 ，M 个 水 听 器 均匀 分 布 为 线性 阵 
列 ， 其 间距 为 A, 满足 人 < 和 /2 空间 采样 条 件 , 其 
中 X = co/f 是 声 源 中 心 频 率 对 应 的 波长 , 9 是 声 源 
的 到 达 阵 列 中 心 的 方向 角 ; 那么 , 声 源 s 的 信号 传 
播 到 水 听 器 m 的 时 间 延 迟 可 以 表示 为 


Asinb0 
Tnm 三 Tn0 十 (m 一 1) 一 = 


(14) 


且 在 [-45°, +45°] 之 间 . 通过 公式 (14) 的 讨论 , 波 
束 形成 的 空间 分 辨 率 仅仅 为 11°, 不 能 分 辨 10° 以 内 
间距 的 两 个 声 源 . 但 是 , Capon 和 MUSIC 方法 都 能 
很 好 地 分 辨 出 两 个 声 源 , 并 且 , Capon 较为 准确 地 分 
辨 出 两 个 声 源 的 功率 , 而 MUSIC 只 获得 了 高 超 空 间 

同样 的 结果 如 图 3 所 示 , 分 别 模拟 了 喷 管 直径 
是 0.1mm 和 5 mm 的 气泡 生成 、 脱 落 情况 , 流速 0.1 
m/s, 气泡 平均 直径 0.3 mm， 网 格 划分 0.1 mmx0.1 
mm。10 个 水 听 器 组 成 2 列 线性 均匀 阵列 , 声 源 成 像 
的 频率 是 1000 Hz， 
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(a) 直径 为 0.1 mm 的 喷 管 气泡 模拟 (b) 直径 为 5 mm 的 喷 管 气泡 模拟 流速 


0 xX/m Xx/m 
.37 0.58 0.59 0.6 0.61 0.62 0.63 
.02 0.57 0.58 0.59 0.6 0.61 0.62 0.63 


医 0.02 
(JJBubble region (NBubble region 
言 Sensor position| 本 Sensor positio 
& 


(9) 喷 管 (0.1 mm) 气 泡 的 波束 (d) 喷 管 (5 mm) 气 泡 的 波束 形成 成 像 


XxX/m 
0.6 


| 0.57 0.58 0.59 0.61 0.62 0.63 
| .02 12 
10 
8 
6 
4 
0 
= 
-4 
| 一 6 
| =8 
(e) 喷 管 (0.1 ao 的 下方 半天 丰 外、 (四 喷 管 (5 mm) 气 泡 的 最 小 方差 无 失真 成 像 
X/m X/m 
0.57 0.58 0.59 0.6 0.61 0.62 0.63 
v0 1 0.0%37 058 059 06 0:61 0.62 0.63 14 
| 0.03 ji 0.03 12 
| 0.04 16 0.04 10 
0.05 14 0.05 
i E i 12 | 有 0.06 
| 全 性 衬 3 
| 一 ”0.07 8 0.07 0 
6 
- 0.08 5 上- 008 -5 
| 0.09 2 0.09| 一 4 
| 0.1 | 0.1 -6 


(g) 喷 管 (0.1 mm) 气 泡 的 子 空间 分 解 的 成 像 人) 喷 管 (5 mm) 气 泡 的 子 空间 分 解 的 成 像 


图 3 不 同 喷 管 直径 下 的 波束 形成 、 最 小 方差 无 失真 、 子 空间 分 解 等 成 像 结果 
Fig. 3 Beamforming, Capon and MUSIC methods for acoustic bubble imaging at nozzles (0.1~5 mm) 
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3 结 论 


本 文 阐述 了 一 种 基于 空间 滤波 器 技术 的 声 源 成 
像 方法 ,不仅 能 从 数量 有 限 的 水 听 器 测量 中 定位 气 
泡 声 源 的 位 置 ， 还 能 够 重 构 声 源 平均 功率 的 分 布 情 
况 。 事实 上 ， 基 于 波束 形成 的 声 源 成 像 是 一 种 稳健 
地 、 快 速 地 和 直接 的 方法 , 它 可 以 被 看 作 是 被 测 声 源 
的 功率 分 布 与 测量 阵列 的 传递 函数 的 卷 积 结果 。 当 
测量 阵列 足够 理想 时 ， 波 束 形成 成 像 即 为 被 测 声 源 
功率 分 布 . 所 以 通常 情况 下 , 波束 形成 的 成 像 分 辩 率 
会 比较 低 、 其 功率 分 布 比 被 测 声 源 的 放大 许多 . 相反 
的 , 子 空间 分 解 方 法 理论 上 可 以 达到 无 限 好 的 空间 
分 辨 率 ， 但 是 其 不 能 重 构 声 源 的 功率 分 布 。 折 囊 来 
看 ， 最 小 方差 无 失真 成 像 方法 既 能 够 提供 足够 高 的 
空间 分 辨 率 、 又 能 够 相对 准确 地 估计 出 声 源 的 平均 
功率 分 布 情况 , 并 且 所 利用 的 水 听 器 阵列 结构 简单 、 
架设 容易 ， 还 能 根据 射流 气泡 的 特性 、 流 速 以 及 喷 
管 形状 和 容器 大 小 来 调整 水 听 器 的 阵列 结构 ， 以 便 
自 适 应 获得 气泡 声 源 的 成 像 结果 . 

在 正在 展开 的 工作 中 , 至 少 三 方面 需 

) 应 用 具有 相关 性 的 单 极 子 或 空间 源 
a 
射流 条 件 ; 

2) 在 密闭 容器 空间 中 , 校 
应 以 便 减 弱 混 响 等 干扰 ; 

3) 针对 在 气泡 生成 、 高 
件 下 的 实际 数据 ,改进 基 于 空 
卷 积 技术 的 声 源 成 像 方 法 . 
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